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Abstract. We propose a general method for finding sharp constants in the imbeddings of the 
Sobolev spaces H"^{M), defined on a n-dimensional Riemann manifold Ai into the space of 
' bounded continuous functions, where m > n/2. The method is based on the analysis of the 

, , asymptotics with respect to the spectral parameter of the Green's function of the elliptic operator 

■ of order 2m, the domain of the square root of which defines the norm of the corresponding 

Sobolev space. The cases of the n-dimensional torus T" and n-dimensional sphere S" are treated 
in detail, as well as some manifolds with boundary. In certain cases when A4 is compact, 
If^ • multiplicative inequalities with remainder terms of various types are obtained. Inequalities 

with correction term for periodic functions imply an improvement for the well-known Carlson 
inequalities. 
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1. Introduction 

In this paper we propose a general method for finding sharp constants in multiplicative 
inequalities of Sobolev-Gagliardo-Nirenberg type characterizing the imbedding of the Hilbert 
Sobolev space H'^{M.) into the space of bounded continuous functions C{Ai). Here M. is an 
n-dimensional manifold and m > n/2. The inequalities are as follows: 

(^ ^\ II l|2 <r T<r\\ 1126* ii ||2(1-^') 

\^-^) W^Wh^iM) - ^\\'^\\h\m)\\'^\\h-^{m)-' 

where -oo < / < n/2 < ?n < oo, so that < 6* = < 1- 

' ' 2(m—l) 

Of course, multiplicative inequalities are known to hold in a much more general Lp setting (see, 
for instance, [3].|25]). when the three norms in (jl.ip are the Lp, Lq and Lj.-norms, 1 < p,q,r < oo, 
and 9 is accordingly defined by scale invariance; in this sense we shall be dealing only with the 
L^o — L2 — L2 case in this paper. In the one-dimensional case, when = M, IR_|_, or (0, L) 
the corresponding interpolation inequalities are called the inequalities for derivatives. There 
exists a vast literature devoted to them, see, for instance, p!].|20j and the references therein. On 
the whole line, the general Lp — Lq — Lj. case of a function and its first-order derivative was 
completely settled in [2Tj . 

The sharp constant in inequality (II. ip on M was found in [23], the more complicated case of 
the half-line M+ was solved in [10] , the value of the constant in closed form was obtained in [15] . 

The sharp constant in inequality (|l.ip for periodic functions with zero mean was found in |12j . 
In particular, for Z = it was shown there that 

ll-ull^ < K(m)||nf^||u('")||2(i-^), n G /7;;j,(0,2^), / u{x)dx = 

Jo 

where m > 1/2, 9=1 — l/(2m) and 



K{m) = (2m9'^{l - 9)^"'^ smTT9^ ^ 



The constant (which is, in fact, the same as in the case of the whole line) is sharp and no 
extremal functions exist. We also observe that the first order inequality, namely, 

(1.2) II^^IlL < Iklllh'lli G ^^pcr(0,27r), / u{x)dx = 0, 

Jo 

was, in fact, proved much earlier in [TT], as a proof of the Carlson inequality. (In the end of this 
section we discuss the connection of our results with the Carlson inequality in greater detail.) 

Sharp constants in inequalities (jl.ip on the sphere §" were found in |13j . 

A comprehensive analysis of sharp constants in (jl.ip and in logarithmic Brezis-Gallouet 
inequalities on the n-dimensional torus T" has been done in [2] , where inequalities with correction 
terms were obtained for the first time. For example, the following inequalities hold 

r-27r 



loo ^ ll^ll II""'!! ll^ll ) u €^ Hp^j.{0,27T), / u{x)dx = 0, 

^ Jo 

r2n 



(1.3) 

||n||2 < ^||uf/2||u''||i/2 - A||^||2 u £ H^^(0,2tt), [ u(x)dx = 0, 

where all constants are sharp and no extremal functions exist. These inequalities have been 
verified in [2] with an essential help of numerics. Purely analytic proof of them is given in the 
present paper as one of the applications, see Theorem 13. 7[ 
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In the first part of this paper we describe a general method for finding sharp constants in 
inequalities of the type (jl.ip including the inequalities with correction terms. Namely, let A 
and B be self-adjoint non-negative elliptic differential operators on Ai of order 2m and 21, 
respectively. (To describe the main ideas we may assume for a moment in this introduction that 
I = and B = Id.) By the elliptic regularity the kernel of A is finite dimensional and consists 
of smooth (orthonormal) functions: 

ker A = span{(^i, . . . , 93^}. 

We set H"^ := H'"^{M) D ker A-^ and := H''{M)n keiA-^. Then the square roots of A and 
B define the equivalent norms on H"^{Ai) and H\JvI): 

\\u\\\ := {Au,u) ~ \\u\\jjrr.(^_M), \\u\\% := {Bu,u) ~ 

for u G H"^ and u ^ H\ respectively. 

We fix an arbitrary point ^ € 7W and consider the following analog of inequality (jl.ip : 

(1.4) \nm' < K\\ua\n\\f-'\ 

In particular, we are interested in the sharp constant K = K{^) in this inequality. To study this 
problem, we introduce the following maximization problem: for every number D > Xq, where 
Ao is the first eigenvalue of B^^/^AB^^/"^ in H := L2{M) n kex , find V(^, D), solving 

(1.5) V(^,I?) := sup||u(OP : ueH"", Iklll = 1, \\u\\a = D 
Then, the inequality 



^111 



holds and, by definition, V is the smallest function for which it holds. Thus, in particular, if the 
maximization problem (jl.5p is solved, finding the best constant K in (|1.4p is reduced to finding 
the smallest K for which the following inequality holds: 

Y{C,D) < KD^-^, D>\o. 

The solution of the maximization problem (|1.5p can be expressed in terms of the Green's function 
of the following elliptic operator of order 2m: 

A{X) := A + \B, A>-Ao. 

Namely, let Gx{x,(,) be the Green's function of it in H: 

k 

A{X)Gx{;0 = ^(-,0, ^(^,0 := S{x,0 -Y.^jix)^,{0, 

where 5(x,^) is the Dirac delta function. From elliptic regularity we see that Gx{-,0 G H"^ C 
C{Ai). Then, as shown in Theorem 12.31 there exist a unique extremal function in (jl.5p 

UD,d^) = TT^ . ^.|| , so that V(e, D) = — 2 ' 

\\Gx{d)[-,^)\\b \\Gx(D){-,^)\\b 

where X{D) is a monotone increasing function with lim£)_>.oo X{D) = 00, lim/)_j.A(, X{D) = — Aq, 
defined as a unique solution of 

\\ud,(\\a = D. 

Moreover, as shown in Theorem 12.51 the sharp constant K in the multiplicative inequality (jl.4p 
can be expressed in terms of the Green's function as follows: 

^ = ^(^) ■= ee^iley-e ■ [>^'Gxic, o} < oo. 
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In addition, the extremal function exists if and only if the supremum with respect to A is 
attained at a finite point A^, < oo and then n*(x) = const Gx^{x,^). We note that the approach 
used in [12], |13j in the case of n > 1 eventually reduces to the same one-dimensional 
maximization problem. 

The further progress (inequalities with correction terms) is based on the knowledge of the 
asymptotic behavior of the Green's function as A — )• oo: 

(1.6) GA(e, = (giX'-' + 92 + 93X'^' + o{X'-')) , 

where gi = 9i{C) > 0, 52 = 92(0 ^ 93 — 93(0 some given numbers. This expansion is 
assumed, and its verification is one of the main technical tasks in particular examples in the 
second part of the paper. 

Given ()1.7p . we prove in Proposition l2.6l that the solution of (jl.Sp has the following asymptotics 
as I? — )• 00: 

(1.7) YiC,D) = 9iSD^'^+g,y-S-^ ^'^^' - ''^^^^ D^-^ + oiD^-^), 3:=^,^^^^. 

Suppose that the third term is negative (this is always the case when 53 < 0), then 

(1.8) V(e,Z))<5i5Z)i-^ + |, 

for large D. We shall see in certain examples in the second part of the paper that this inequality 
holds for all D. This implies the multiplicative inequality with correction term (which is negative 
if(72 <0): 

I lc\\2 ^ Q\\ \\Q II ||l-6» I 92 II ||2 

\u{i)\ <gxS\u\j^\u\j^ +-q\Ab 

with best possible constants. 

In the end of first part of the paper we give in Theorem 12.81 the following variational charac- 
terization of V(^, D): 

V(^,D)= iiif {(A + Z5)Ga(^,0}. 

Ae[— Ao,oo) 

This formula is very useful for proving in certain cases (which are few!) that (jl.Sp holds for 
all D. The scheme is as follows. Usually it is impossible to find the unique minimizer A(D) 
explicitly, but we always have the asymptotic formula for it 

A(D) = -^D + ^51)^ + . . . , 

see (|2.44p . We somehow truncate this expansion and denote the result as \^(D^. Then to 
prove (jl.Sp we proceed as follows 

N[i,D)-g^SD^-' - I < (A.(Z)) +Z))G,,(B)(e,0 - |. 

Now for a fixed the right-hand side is an explicit function of D only and more or less standard 
estimates can be used prove that it is negative for all D. In fact, inequalities (jl.Sp as well as a 
number of other inequalities with lower order correctors mentioned below are verified using this 
scheme. 

In the second part of the paper we consider examples and applications of the general approach 
described in the first part and we use the above scheme in Theorem 13.71 for purely analytic proof 
of inequalities (jl.Sp (the proof in [2] involves some reliable computer calculations) . 

We first deal with manifolds without boundary and consider the case A = (—A)™, B = (—A)', 
where A is the Laplace-Beltrami operator. Then H^{M) = H^{M) n {/^ u{x)dM = 0}. The 
asymptotic formula ()1.6p for the Green's function on the torus is obtained by the Poisson 
summation formula. This works when 1 = 0. For / > there is a singularity, which can be 
removed by differentiation. This, in turn, produces the problem of finding integration constants. 
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The corresponding technique was proposed in and it is further developed here. We consider 
only one example n = 3, m = 2, / = 1 and prove the following inequality on the three-dimensional 
torus T^: 

\\u\\lo < 7?-||V7x||||A'u|| - ^\\Vuf, u G H'^{T'^), 

where the integration constant /Ja is expressed in terms of a super-exponentially convergent 
series, (3^ = —8.91363291758515127. Both constants are sharp and no extremals exist. 

Next we study inequalities on spheres. On ^e consider the case when A = (—A)™, m > 1/2, 
and B = I. The corresponding Green's function is independent of ^ and is given by the series 

1 °° 1 
G'(A) = — /i"^5^(2n + l)<^Mn + l)), where ^i = \-^'"', '^^'^^ = 

n=l 

Thus, we need to find the asymptotic behavior of functions of the type 

oo 

F(^) = ^(2n + l)/(^n(n + 1)) as /U ^ 0, 

n=l 

where / is sufficiently smooth and sufficiently fast decays at infinity. This is achieved with the 
help of the Euler-Maclaurin formula in Lemma 13.141 

I 2 1 

F(^) = - / f{x)dx - -f{0) - -^/'(O) + Oifi'). 

II Jq 3 15 



This gives the asymptotic expansion of the type (jl.6p for the Green's function and, hence, the 
asymptotic expansion (|1.7p for the solution Y{D) of the corresponding maximization problem, 
in which the second and the third terms turn out to be both negative. Therefore a negative 
correction term may exist. For m = 2 we show that this is indeed the case and the following 
inequality (with two sharp constants and no extremal functions) holds for u G H'^{S'^): 

„ „2 1 , „ 1 „ „2 

For larger m [m = 3, m = 4) the negative correction terms still exist, but are smaller than the 
second terms in the expansion for Y{D). 

On we consider only one example with A = (— A)^ and B = —A. The series expressing 
the Green's function can be summed in closed form, and along the same lines we obtain a sharp 
inequality for u G H'^(E>^): 

\\u\\l<l.\\Vu\\\\Au\\-^jVuf. 

In the remaining part of the paper we consider manifolds with boundary. We first prove a 
sharp multiplicative inequality on the half- line for the Bessel operator |17| . 

Then we consider the case when u G Hq{0,L). The correction term still exists, but is expo- 
nentially small, namely, the following inequality holds: 

Iklloo < ll^^lllk'll(l - 2e ). 

Both coefhcients on the right-hand side are sharp and no extremal functions exist. 
For a second order inequality on the interval (0, L) 

II l|2 ^ T<'\\ 11^/2 II //||l/2 
ll-Ulloo ^ -'^ II^IIl2(0,L)II'" IIl2(0,L)' 

in the case when u G Hq{0,L) we can use extension by zero, and therefore K is the same 
as on M, namely, K = (see, [23], and also (j3.2p ). In going over to a wider space u G 
H^{0,L) n H^{0,L) the constant may increase. In Theorem 13.211 we show that this is indeed 
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the case: K = 4^ • coth ^ = 4^ • 1.09033 . . . , and, in addition, there exists a unique extremal 



function. 

A somewhat opposite result is obtained in Theorem 13.231 the constant ^ in the inequality 

\\u\\l^ < ^||Vu||||An||, 

is sharp both for u G H^{n) and u G H^{n) D H'^{n), n C M.^ . 

In conclusion we observe that inequalities of the type (jl.ip have applications to various prob- 
lems in partial differential equations and mathematical physics. For example, the far-going 
generalization to the matrix- valued case [6] of the simplest inequality 

< hlllk'H, ueH^{M) 

gives the best-known estimates of the Lieb-Thirring constants for the negative trace of the 
Schrodinger operators in M" [19| . 

Accordingly, inequalities with correction terms (II. 3p imply by the method of [7] a simulta- 
neous bound for the negative trace and the number of negative eigenvalues for the Schrodinger 
operators on [14] . 

Finally, inequalities (II. 3p provide an improvement to the well-known Carlson inequality [5] 
(see also |18] and the references therein): for Ofc > 



(1.9) i^a, <WE«n 



\k=l / \A:=1 / \fc=l 



where the constant vr is sharp and the strict inequality holds unless all = 0. In fact, as shown 
in [11], inequalities (jl.2p and ()1.9p are equivalent, hence an improvement to ()1.2p results in an 
improvement to ()1.9p . 

Given a sequence > 0, we set ao = and consider a real-valued function 



with mean value zero, for which 

oo oo cc 

ll^^lloo = ^^(0) = 2^0fc, lliilp = 47r^ a|, = 47r ^ fe^o; 

fc=l k=l k=l 



Now inequalities (ll.3p can be reformulated as follows 

(oo \ 2 / oo \ 1/2 / oo \ 1/2 oo 

^ ^ k=i / \k=i / \k=i / k=i 

/ oo \ 2 /— / oo \ 3/4 / oo \ 1/4 „ oo 



\fc=l / ^ \fc=l / \fc=l / k=l 

where all constants are sharp and the inequalities are strict unless all = 0. Moreover, the 
interpolation inequalities on the tori T" discussed above can be naturally considered as multi- 
dimensional generalizations of the Carlson inequality. 

In this paper we use standard notation. Sometime the C-norm (or Loo-norm) is denoted by 
II • I loo and the L2-norm is denoted by || • ||. 
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2. Part I. General theory 

2.1. Assumptions and preliminaries. Let M be an n-dimensional compact Riemann mani- 
fold and let 

{u,v) := / u{x)v{x) mes{dx) 

Jm 

be the standard scalar product in H := L^(A^,mes) (where mes stands for the measure on 
A4 associated with the Riemann metric). As usual, we define the Sobolev spaces W^'^{A4), 
1 < p < oo as spaces of distributions whose derivatives up to order I belong to L^^Ai) (this 
definition works for integer I only, for non-integer and/or negative I, the spaces W'''^{Ai) are 
defined in a standard way using the interpolation and duality methods, see e.g., [25] for the 
details). For the case p = 2, we will denote by H"^ = H^{M) the Sobolev space W^''^{M) 

We assume that A is an elliptic self-adjoint differential operator of order 2m on M with 
smooth coefficients. This operator is supposed to be non-negative 

(2.1) (Au,u)>0, ueH"^, 

although it may have a non-zero kernel. Then, due to the classical elliptic theory (see e.g., [25j). 
the kernel is finite-dimensional and is generated by smooth functions 

(2.2) ker^ :=span{99i,--- (^i G C°°(A1), Aipi = 0. 
Without loss of generality, we may assume that 

||y?i||/f = l and {ipi,ipj) = 0, i^j. 
We denote by H the orthogonal complement of ker A in H = {Ai ) and define 

(2.3) := H'{M)nH, s G R. 

Then, due to the elliptic theory, A is an isomorphism between ^■^+2'" ^nd for all s G M and, 
in particular, the equivalent norm in H"^ is given by 

(2.4) \\u\\\:= {Au,u), ueH"^. 

We also introduce the second elliptic non-negative and self-adjoint differential operator B of 
order 21 < 2m on Ai with smooth coefficients such that 

(2.5) kerScker^, and S ker Acker A. 

Then, as not difficult to see, the operator B is an isomorphism between ^^+2' and H"^ for all 
s G M and, in particular, the equivalent norm in is given by 

(2.6) \\u\\l := {Bu,u), ueH^. 
In addition, we have the analogue of the Poincare inequality: 

(2.7) IKIll < Aoinll^, uG^™, 
where Aq > is the minimal eigenvalue of B~^/'^AB~^/'^ in H. 

2.2. Interpolation inequality and associated variational problem. If 

Tl 

(2.8) I < - <m, 

then for every u G H"^ the following interpolation inequality holds: 

(2.9) \Mc(M)<C\\u\\U\u\\]^', ^ = W + 7n{l-e), 
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see [3], [25]. Our aim is to refine inequality (j2.9p and, in particular, to find the best constant 
C = C(n, m, or/and the lower order extra terms in it, etc. To this end, we fix an arbitrary 

point £ Ai and a positive D and consider the following maximization problem: 



(2.10) V(^,D) := sup||u(Or : ueH"", I^IIb = 1, 1^111 = ^ 

Indeed, in view of (12. 9p and (j2.7p the function V is well-defined for aW ^ £ M and all D > Xo. 
On the other hand, due to the homogeneity. 



(2.11) |^(^)|2<||^|||v(^^,Q^j, ee-M, ueH^. 

In particular, the best constant C in ()2.9p is the minimal one for which the inequality 

(2.12) supV(C,-D) < C^D^"^ 



|2 



holds for all D > Xq. This reduces the study of inequality (12. 9p to the investigation of the 
maximization problem (j2.10p . 

2.3. Green's functions and reproducing functionals. In this subsection we prepare some 
technical tools which are necessary to give the analytic description of the function V in terms 
of the Green's functions of the appropriate elliptic operators on M. Namely, let 

(2.13) A{X):=A + XB, A > -Aq, 
and let the function Gx{x,(^) solve 

k 

(2.14) A{X)Gx{;0 = K,0, 5(x,e) = 5(x,0-E<^^(^)'^*(^)' 

1=1 

where 6{x,(^) is the Dirac ^-function at point ^ G and 6 is its 'projection' on the space H. 
Then, due to the Sobolev embedding theorem and (|2.8p . 

^(•,6g^"'" 

and, therefore, due to the elliptic regularity, 

(2.15) GA(-,Oe^"cC(A^). 

In addition, since A(A) is self-adjoint, Gx{x,^) = Gx{^,x) and, consequently, (12.15j) implies that 

Gx G GiM xM). 

Moreover, 

(2.16) uiO = (^(., e), u) = iA{X)Gx{;0,u) = (Ga(-, 6, A(A)tx). 
for all u £ H^. In particular, taking u{x) := Gx{x,£,), we have 

(2.17) GxiCO = {A{X)Gxi;0,Gx{;0) = I|Ga(-, OIIa + A||Ga(-, Olll > 0. 

The following simple lemma is nevertheless the main technical tool for the method of reproducing 
functionals and will allow us to find the analytic expression for the function V. 

Lemma 2.1. Let the above assumptions hold. Then, for every u G and for every A > — Aq 

(2.18) < GxiCO {\Ma + MMl) ■ 
Moreover, the equality holds if and only if u{x) = cGx{x,C} fof some c G M. 
Proof. Indeed, since A(A) is positive definite, by the Cauchy-Schwartz inequality, 

lOP = (Ga(-,0, A(A)n)2 < (A(A)GA,GA)(A(A)n,n) = Gx{tO {Ml + MMl) 
and the equality here holds if and only if u{x) = cGxix,^,). Thus, the lemma is proved. □ 
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We note that, up to the moment, A > — Aq is a free parameter in (j2.18p . The next lemma 
shows that this parameter can be chosen in such way that the quotient ||Ga(-)|Ia/II^a(") Oils 
achieves any prescribed value. 

Lemma 2.2. The function 
(2.19) 

is strictly increasing on [— Ao,oo). Moreover, 

(2.20) lim D{X) = Aq and lim D{\) = oo, 

A— > — Ao A— >+oo 
SO the inverse function D — t- X{D) is well defined on [Ao,oo). 

Proof. We will show that D'{X) > 0. To this end, we note that the function G\{x,S,) is smooth 
with respect to A and the derivative G';^(x,^) := dxGx{x,^) solves the equation 

(2.21) AG'x{;0+>^BG'x{;0=H>^)G'x{;0 = -BGx{;C), 
so that 

AG'xi;C) = -BGxi;C)-XBG'x{;0, G'x{;0 = -A{X)-'BGx. 

Therefore, 

(2.22) l\\G,{;0\\%D'{X) = 

= \\G,{-,mUAG',{-,C),G,{-,0)-{BG',{-,C),G,{-,mGx{-,m\ = 
= -\\Gx{;m% - X\\G,{;0\\l{BG',{;0,G,{;0) - {BG',{;0,G,{;0)\\Gx{;0\\a = 
= -\\Gx{;0\\% - {BG',{;0,Gxi;OmX)Gx{;0,Gx{;0) = 

= {A-\X)BGx{;0,BGx{;mH^)Gx{;0,Gx{;0) " I|Ga(-,OI|Ij- 
By the Cauchy-Schwartz inequality we have 

(2.23) \\v\\% = {Bv, vf = {A-'^/\X)Bv, A^/\X)vf < 

< \\A^^/^{X)BvfH\\A^^\X)v\\jj = {A^H^)Bv,Bv){A{X)v,v) 

for all V G H^. Taking v{x) = Gx{x,£,), we see from ()2.22p that D'{X) > 0. Moreover, equality 
in (|2:23]1 holds only if 

aA{X)v = f3Bv, q,/3gR 

and from elliptic regularity v should be at least C°°-smooth. However, v = Gx{-,0 cannot be 
C°°-smooth, so the equality is impossible and D'{X) > 0. Thus, we have proved that D{X) is 
strictly increasing. 

We now need to verify (|2.20p . We start with the limit A — )■ — Aq- To solve (j2.14p near A = — Aq, 
we introduce w{x, ^) = B^/^Gx{x,^). Then 

(2.24) B-^/^AB-^/^w + Xw = B-^'H{x, i) 

and, by definition, Aq is the smallest eigenvalue of the self-adjoint positive operator B~^/'^ AB~^/'^ 
in H. Clearly, this operator has compact inverse, so its spectrum is discrete. Let {V'fcjfcLi be 
family of the orthonormal eigenfunctions of it, and let the first p {p > 1) of them correspond to 
the smallest eigenvalue Aq. Clearly, ipi £ C°°(M). Seeking w in the form 

p oo 

(2.25) w{x,^) = '^ai7pi{x) + ^ aiipi{x) =: wq + , 

1=1 i=p+l 
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substituting this into (j2.24p . and taking the scalar product with ipi, . . . jipp, we find the coeffi- 
cients of the WQ-pait of the solution: 



A + Ao^ V-.'-. ^^^^ 

where we set := B~^/'^'il)i. Observe that the t(;"'"-part of the solution remains bounded as 
A ^ -Aq. Applying B~^l'^ to H^TM we find Gx{xX)-- 

1 ^ 

(2.26) Gx{x,i) = -— -^^,(0^,(x) + Gf(x,0, 

A + Ao^ 

where the part Gj^{x,£^) remains bounded as A — )■ — Aq- Thus, 

IIELi*«)*(OIII Ef.i*({)^ 

Let us consider the case A — )■ oo. Assume that ()2.20p is wrong and we have 

lim D{\) = Dmax < oo. 

A— >oo 

Then, 

(2.27) \\Gx{■,0\\A<D„^ax\\Gx{■,0\\l■ 

Multiplying equation (I2.14p by G\{x,S^) integrating in x G , and using the embedding ^™ C C 
and l^in . we get 

\\Gx{;mA + M\Gx{;0\\l = GxiC,C) < I|G'a(-, Olb < C||GA(-,e)IU < Ci||Ga(-, Ob- 
Due to this estimate, we have 

\\Gx{;mA<C\\Gx{;0\\B<C2X~\ 

where C2 is independent of A — 00. Thus, wxix,^) := AGa(x,^) is uniformly bounded in H"^ 
and, without loss of generality, we may assume that wxi-jS,) — ?• i'^oo(")'^) weakly in this space. 
Then, obviously, w 00 € H"^ and 

In particular, since Woo G H"^, we have € Due to the assumption (j2.8p we have 

m — 21 > —n/2, so that (^(•,i^) € for some s < n/2 which is impossible. Thus, (|2.20p is 
proved and the lemma is also proved. □ 

2.4. Main result. The aim of this subsection is to give the analytic expression for the function 
V in terms of the Green's functions introduced above. This result is stated in the following 
theorem. 

Theorem 2.3. Let the above assumptions hold. Then, for every D G [Aq, 00) and every ^ ^ Ai, 
the supremum in (j2.10p is the maximum and this maximum is achieved in a unique point 

(2-28) UD,^{x) := — , 

\\Gx(d)[-,^)\\b 

where the function X{D) is defined in Lemma \2.SX In particular, 

Gx{D){LO 



(2.29) Y{C,D) 



\GxiD)i;0\\ 



B 
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Proof. Indeed, let u E H"^, be such that \\u\\b = 1 and = D. Then, according to Lemma 
O there is a unique A = X{D) which solves (j2J9|) . According to (fZTS]) with A = X{D), 

<Gx(D){tO{D + XiD)) 

and the equality here holds if and only if u{x) = cGx[d) i^y 0- Taking the i?-norm of both sides 
of this equality, we see that c = ||G';^(d)("5 0IIb^ u{x) = ud,^{x). This finishes the proof of 
the theorem. □ 

Thus, in order to find V, we need three functions 

(2.30) /^(A) := GA(e,e), 55(A) := \\Gx{;0\\l, ■= \\Gx{;0\\a- 
Then 

(2.31) D{X) = and V(e, A) := V(e, Z)(A)) = X G [-Aq, 00) 

5? (A) 5? (A) 

and we have the parametric representation of the function D — )• V(^, D). The next lemma shows 
that the functions g and h can be expressed in terms of /. 

Lemma 2.4. The functions f, g and h satisfy the following equalities 

(2.32) /^(A) := i-/^(A) = -5^(A), /15(A) = /g(A) + A4(A). 

Proof. Indeed, multiplying (j2.14|) by G\{x,^) and integrating it over x G A^, we have 

(2.33) /^(A) = GA(^,e) = I|Ga(-,6IIa + A||Ga(-,OII| = + A5s(A). 
Differentiating this formula with respect to A and using (|2.2ip . we get 

(2.34) 4(A) = 2(A(A)G';,(-,e),G'A(-,0) +5?(A) = 

= -2(i3GA(-,0,GA(-,0) + 5C(A) = -2<7c(A) +5c(A) = -<7c(A). 
It only remains to note that (j2.34p and (|2.33p imply ()2.32p and finish the proof of the lemma. □ 
The next result shows that the sharp constant in the inequality 

(2.35) |n(OP < K\\uf^\\uf}^~^\ u G if™. 

is expressed in terms of the following scalar maximization problem again involving the Green's 
function G\{£,.,^). 

Theorem 2.5. For a fixed ^ G the constant K in (j2.35p is given by 

(2.36) K = KiO := ge^ilgy^e ' f^^ [>^'Gx{i, o} < oo, 

where 9 is defined in (j2.9p and the constant (|2.36p is sharp. Furthermore, the extremal function 
in (j2.35p exists if and only if the supremum in (j2.36p is attained at a finite point A* . 

Proof. We first show that ()2.36p is finite. Indeed, using ()2.17p and the inequality (12. 9p (with the 
non-optimal constant C!), we have 

X'G.iCCf < X'\\G,{C,-)\\loo < C{X\\G,{C,-)\\l)'{\\Gx{C,-)fA)'~' < 

<c(iiGA(e,-)iii + A||GA(c,-)iil) = cGA(e,e) 

and ()2.36p is finite. 
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Let us check that (j2.35p holds with K = K((^). Indeed, let u S H"^ be arbitrary and let 

II II 2 

A := . Then, using (j2.18p . we have 

(2.37) HOI' < Gx{C,OMl (p| + ^) = ^GA(e,OA||n||| = x{t >^'^ X'^' \\u\\l = 

Let us check that K{^) is sharp. We first assume that the supremum in (|2.36p is the maximum 
which is achieved at A = A*. Then, 

= ^{\'h{mx=K = At(4(A.) +^A;Vc(A*)) 
and A/^(A*) + 61 /^( A*) = 0. Therefore, due to (l2:32D . 

^^(^*) _ /g(A0 + A.4(A.) _ i-e ^ 

f) \\G (£ )\\'^ 

Thus, A=^ = j^ uril'-w^ inequalities in (fOTl) become equalities if we take u{x) = 

Gx,{^,x), SO we have the exact extremal function in that case. 

Since X^G\{^,^) — )• as A — ?• 0, we only need to consider the case when the supremum in 
()2.36p is achieved as A — t- oo. Then, two alternative cases are possible: 

1) there are sequence {A^j^j^ of local maximums such that 

AfcGA,(e,e)^sup{A^GA(e,e)} 

A>0 

Since the derivative vanishes at local maximums, then arguing as before, we see that the sequence 
of conditional extremals := Gx^{(,,x) does not allow us to take the constant K strictly 

less than K{^) and (12.36P is sharp. 

2) The function A^Ga(C)C) is eventually monotone increasing as A — t- oo. Then the limit 

(2.38) Goo := lim A^GA(e,e) 

A— >oo 

exists and is strictly positive. Using the fact that the derivative is integrable we can find 
sequences Afc — )• oo and — >• such that 

|^(AVc(A))|A=A,=e.A,T^ 

This, together with ()2.38p gives 

A^/^(Afc) = Goo + OA^oo(l), A^+V^(Afc) = -0Goo + OA^oo(l). 

Therefore, 

(2.39) — = — h OA->oo(l)- 

Finally, taking Uk{x) := Gxi^{£,,x), after straightforward transformations we see that 

(2.40) ., .,b«>lL., ^teJ.f^V-°.fl.°'^>i 



Passing to the limit A; — t- oo and using (j2.39p . we obtain exactly K{^) in the right-hand side and 
verify that -fC(^) is sharp in the second case as well. 
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To complete the proof it remains to show that if there exists an extremal function in (j2.35p , 
(j2.36p . then the supremum with respect to A in (|2.36p is attained at a finite point. 
Using the elementary identity for positive a, b 

Xe^2ii-e)f^20 ^ ^.(1 _ ^)(i-.) (^2 + ^^^2) ^ _±_^^ 

we have for the extremal function the equality 



= supl X'^Gx{^,0\Ta {\W*\\a + ^*\W*\\b) ^ where A* = — 
A>o I J 1 

Lemma |2. II now gives that necessarily n*(x) = Gx^{x,S^) and 



1^* II A 



1^* 111? 



sup|A^GA(e,e)|i = GA.(e,0, or sup (A'^GAie, o| = A^Ga. (C, 0- 

A>0 I ) K A>0 I J 

The proof is complete. □ 

2.5. Asymptotic expansions for big A. In this subsection, we derive some useful formulas 
for the function V(^, D) when D — )• oo. To this end, we need to know the asymptotic behavior of 
the Green's function /^(A) = G\{^,^). It is not difficult to see using the localization and frozen 
coefficients technique that the limit (j2.38p exists and is strictly positive, so the leading term 
in the asymptotic expansions of Gx is known. However, the further terms in the asymptotic 
expansion seems problem dependent and we do not know the general formulas for them. By this 
reason, we just assume that 

(2.41) /^(A) = A-i (51 A^-^ + 52 + <73A^-' + o{X'~') 

where ffi = 51(C) > 92 =52(0) 53 = dsiO some given numbers, and that we are able 
to differentiate the expansions (I2.4ip with respect to A. This assumption will be satisfied in 
most part of our applications. Then, the following result holds. 

Proposition 2.6. Let the above assumptions hold and let, in addition, the asymptotic expansions 
()2.4ip be true. Then the following asymptotic expansion holds as D ^ 00: 

(2.42) Y{C,D) = giSIji-^+g2^-lg-i ^2(l - 2 - ^^gig3 ^e-i^^(^g-i)^ ^ 1 



2 03^1 ^ e(>{i-ey~s- 

Proof. The proof of this proposition is a straightforward (although rather technical) computa- 
tion. Setting for brevity gi = a, g2 = b, g^ = c we have 

/(A) = aX~' + 6A-1 + cA-(2-^) + o(A-(2-^)), 

5(A) = a0A-(i+^) + 6A-2 + c(2 - 0)A-(3-^) + o(A-(3-^)), 

h{X) = /(A) - Ag(A) = a(l - e)X-^ - c(l - 0)A-(2-«) + o{X'^^~'^'^). 

Next, we find the asymptotics, as D — )• 00, of the (unique) solution A = X{D) of the first equation 
in (|23T]1 : 

h{X) A(a(l - 6*) - c(l - 0)A-(2-2e) + . . . ) 



D{X) 



g{X) aO + 6A-(i-^) + c(2 - 0)A-(2-2^?) + . . . 
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or 

(2.43) X-AX<^ + BX-'+'' + ... = 5:=D-^^, A = A, B=^^^^. 

I — a ad a^y^ 

The unique large solution of this equation has the asymptotics as 5 — )• cx) 
where we find C by substituting the last expression into (j2.43p . which gives 

or, finally, 

(2.44) X{D) =rD + sD^ + tD"^^'^ + . . . , 



where 



^_b 1 ^ _ 2ac9 - (1 - e)b'^ 



It remains to substitute H^M^i into V = p/g = D ■ p/h, for which we have the expansion 

For each power we obtain from p.44p . respectively, 

w N 1 f) ^0(6 + l)s2 - 29tr n , 

A(Z))-i = Vi--i 1)^-2 + ..., 

A(Z?)-^ = ^I)-^ + .... 

Substituting this into (j2.45p . multiplying by D, we obtain after quite a few miraculous cancel- 
lations the asymptotic expansion ()2.42p . □ 

Remark 2.7. We see that, in particular, if for some M. the third term in (j2.42p is negative 
(this is always the case when 53 < 0), then 

(2.46) V(e,I?)<r7i5i^i-^ + |, 

for large D. Then, using the numerics which is reliable for relatively small -D, we will show that, 
in some cases, inequality ()2.46p holds for all values of D. This will give us the improved version 
of D: 



(2.47) \u{^)\^ <g^S\\ufM^^' + %\\nfB 

with best possible constants. 

We also mention an interesting fact that 



D'(X) 



7g(A)+A/^(A) y ^ 2{r^(x)f-ri(x)f(x) 



and, therefore, the positivity of D'{X) proved in Lemma 12.21 is equivalent to the strict convexity 
of the function A — )• j^j^- 
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2.6. Variational characterisation of V(^, A). In this subsection, we give a simple, but very 
useful description of V(^, A) in terms of the Green's function Gx{(,, which does not involve the 
derivatives in A and which allows us in many cases to prove inequality (j2.46p analytically for all 
admissible values of D. Namely, the following theorem holds. 

Theorem 2.8. Let the assumptions of Theorem \2.3\ hold. Then, the function V(^, D) defined 
in (|2.29p can be expressed as follows: 

(2.48) n^,D)= inf {{X + D)Gx{^,0}, 

Ae[— Ao,oo) 

where D G [Ao,oo) and the Green's function Gx{C,0 ^-^ defined in (|2.14p . 

Proof. We first note that, due to the results obtained above the function A — )■ (A + D)Gx{$,,(,) 
tends to +oo when A — t- oo or A — ?• — Aq for every D € (Ao,oo). Thus, the infimum in ()2.48p is 
achieved at some point Aq = Ao(-D) inside the interval. Obviously, this Aq solves the equation 

see (j2.30p . ()2.3ip and (j2.32p . Thus, the equation on the minimal value Ao(-D) coincides with 
equation (j2.19p for A(-D) which has a unique solution and, therefore, the minimum in (j2.48p 
is achieved exactly at Aq = A(-D). It remains to check that the value at this point is exactly 
V(.^, A). To this end, we note that at the extremal point the numbers D and A satisfy the first 
equation of (j2.3ip and, at that point 

9d^) 5? (A) 

where we have used (j2.3ip . Theorem 12.81 is proved. □ 

2.7. Generalizations. Here we briefly discuss several possibilities to relax the assumptions on 
the manifold A4 and operators A and B. 

I) The operators A and B should not necessarily be elliptic differential operators. All the 
theory works word for word if we assume that A and B are elliptic pseudo-differential operators 
(for instance, A = (—A)™ and B = (—A)' where A is the Laplace-Beltrami operator on A4) and, 
in particular, the numbers m and / may be not integers. This allow us to study the interpolation 
inequalities in fractional Sobolev spaces. 

II) The theory can be naturally extended to the case where the manifold Ai has a boundary. 
In this case, we need to assume that the elliptic operators A and B are endowed by the proper 
boundary conditions and these boundary conditions are chosen in such way that 

(2.49) V{A) C V{B), 

where T>{A) and T>{B) are the domains of the operators A and B. It is not difficult to see 
that under this extra assumption(s) the above developed theory remains true for manifolds with 
boundary as well. 

III) The spaces may be further restricted, for instance, we may consider not all functions, 
say, in a disk, but only radially symmetric ones. If the operators A and B are also radially 
symmetric, all the theory works in this case as well. 

IV) The above theory works in many cases where the manifold M is not compact. The only 
problem here is that, unlike the compact case, the Green's functions Gx{x,£,) may have bad 
behavior as x, ^ — )• oo and as a result, the integrals used above may not have sense. So, in 
general theory, one should be accurate with the extra assumptions on the non-compact part 
of M and operators A and B as well as with the possible continuous spectrum of A and B. 
However, in all our "non-compact" applications, these questions will be obvious and transparent. 
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SO in order to avoid the technicalities, we do not present here any "general theory" for the non- 
compact case. 

3. Part II. Examples and applications 

3.1. The case of M". To illustrate our method we consider the simplest case when M = M", 
and let / and m satisfy — oo < I < n/2 < m < oo, so that < 6 = 2(11^) ^ 

Theorem 3.1. The following inequality holds 

(3.1) II^IIL < CM.(Z,m)||(-A)V2n||2e||(-A)-/2^f(i-'^), 

where the sharp constant c^n(l^m) is 

a{n)TT 



(3.2) CR.(/,?n) - (27r)- j^^ ^^^^21 ^ _ 0^\^\2m " (27r)-2{m - l)eO{l - 9)^-0 ^mTrO ' 
Proof. Using the Fourier transform /(^) = J-'(/)(C) = (27r)^"/^ J^„ e"^^^ f{x)dx, for A > 

NIL < (2vr)-" f / mmy< 

(3.3) (27r)-"/ |S(O|'(A0|^|'' + (l-e)|eP 



1 2m 



d^ 



di 



A0iep' + (i-0)iep"^ 



Next, setting A = A* = || (-A)"^/^^^ /II (- A)'/2n||2 we have 

(3.4) Ke\\{-A)^l^uf + (1 - 0)||(-A)™/2^f = A^||(-A)'/2uf^||(-A)™/2^f (1"^), 

which gives 



A 



d^ 



||(-A)'/2^x||2^||(-A)™/2^x||2(l' 



A.^|e|2' + (l-e)|^|2™_ 

The expression in brackets is, in fact, independent of A*, which gives (j3.2p . The fact that the 
constant is sharp can be verified by substituting 

(3.5) ux{x) = F-^ U\ei\i\^' + e2\i\^^'-^^ 



and calculating the corresponding integrals. A simpler way, however, is to observe first that for 
u\{x) all inequalities in ()3.3p become equalities. Next, differentiating the expression in brackets 
with respect to A, one can see that A = ||(— A)'/2ii;^|p/||(— A)'"/2ii;^|p, which proves that (|3.1|) 
for u = u\ becomes an equality. □ 

Remark 3.2. Of course, we get the same result by applying Theorem 12.31 Here 
and 

r (F F\ = \-s 

^AV?,?; A (2^)„2(m-Osin^0 ■ 
Corollary 3.3. For any D > the maximization problem 



(3.6) ¥(1?) := sup<^ |n 



||n(0)p: \\{-AfM' = l, ||(-A)-N|2 = ^| 



has the solution 

(3.7) V(D) = CKn(/,m)i)i' 



GREEN'S FUNCTION ASYMPTOTICS AND INTERPOLATION INEQUALITIES 17 
The unique extremal function is 



see 

Remark 3.4. The integral was calculated using the formula 

oo 

t'ifi\ f ^ 1 p f m+l m+l 



which will also be helpful in what follows. 

Remark 3.5. In the ID case the constant ciR(/,m) was found in |23j . 

3.2. Symmetric manifolds. In this section, we discuss the case where the underlying manifold 
JH is symmetric and the operators A and B are invariant with respect to the symmetry group. 
As a consequence of this, the Green's function Gx{x,^) introduced in (j2.14p depends only on 
X — (since every two points on Ai can be identified by the proper symmetry map). By this 
reason, the key function Gx{(,,(,) is in fact independent of and is a function of one variable A: 

(3.9) Gx{tO = GiX). 

This observation simplifies greatly the analysis and allows us in many cases to compute explicitly 
the best constants in the appropriate interpolation inequalities. We restrict ourselves below 
only to consider the two model examples: tori and spheres although the developed technique is 
applicable to other symmetric manifolds as well. 

3.2.1. The tori. We start with the case of n-dimensional torus T" := [— vr, vr]" and the inequalities 
of the form 

(3.10) II^IIl^(T") < Q,„„„||(-A)'/2n||i^(Tr„)||(-A)-/2n||i;(V) 

for the periodic functions u S if"*(T") with zero mean ((u) := Jj„ u{x) dx = 0). In that case, 
A := (—A)"*, B = (—A)', ker^ = ker B = {const} and all of the assumptions of the above 
developed abstract theory are satisfied if / < ^ < m. In particular, 

(3.11) H"' = //"(¥") n{{u) = 0}. 

The Green's function of the operator A(A) = (— A)™ + A(— A)' on the torus T" (with zero mean) 
can be found by expanding it into the multi-dimensional Fourier series. This gives 

1 gjfc(3;-?) 1 1 

(3.12) Gx{x,0 = -^X^^J^^;^^^J^, GA(e, = G{X) = ^ ^ , 

where k = {ki, • • • , kn) is the multi-index, \k\'^ = kf + ■ ■ ■ + k'^ and the summation holds for all 
multi-indexes /c = (fei, • • • , fe^) / (0, • • • , 0): ZJJ = Z" \ {0}. 

3.2.1.1. The case / = 0. The particular case I = has been studied in [2]. In that case, the 
asymptotic behavior of G(A) as A — )• oo can be analyzed in a straightforward way using the 
Poisson summation formula (see, e.g., [22]): 

(3.13) Yl /W^) = (2^)"^V" E /(2^"^/^)' 

where J^(/)(0 = fiO — (2vr)~"/^ J^„ f{x)e~^^'^dx is the Fourier transform and fi > 0. 
Namely, let 

$a(^) = (2^)-"/2-F-^(l/(|eP'" + A)) = —V f i5o ^"^1 
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be the fundamental solution of the differential operator A(A) in the whole space M^. Then (j3.13|) 
gives 

(3.14) G(A) + — i— A^^ = — i— y -p^^^ = y $A(2vrA:). 

Furthermore, due to scaling invariance, ^x{x) = A~"'^A'"/''^''"^<&i(A^/P''"-*x) and, due to the ana- 
lyticity of the function ^ — )• np^-jiY, we have |<^i(rE)| < C'e"*-^"'^' for some positive C and C^. 
Therefore, (j3.14p reads 

(3.15) G(A) = A'U"/(2m)^^(Q) _ ^^A^i + 0(6-^-^'^^'"'). 

(27r)" 

Finally, computing the table integral '&i(0) = (27r)~" Jjg„ j^^-^, we end up with 
,3,16) G,A) . (2.)-"A- (^^^A* - l) + O,.— 

where cr{n) = 27r"'/^/r(?i/2) is the surface area of the (n— l)-dimensional unit sphere. Note that 
1 — 9 = n/2m now, so (I3.16P has the form of (j2.4ip with (73 = and, due to 12.61 we have 

(3.17) Y{D) = c„(m)D"/(2"') - fe„(m) - /^(m)^-"/^^™) ^ o(D-"/"^), 

where 

(3.18) 

Tra{n) 2m 
cn{m) := CMn(0,m) = (2vr)«nW(2m)(2^ _ n)i-V(2m) - ' ^"^"'^ "= (27r)"(2m - n) ' 

2„l+n/(2m) 2 i |n 



(27r)"7rc7(n)(2n - m)2+"/(2m) ' 
see [2] for the details. Thus, since /n("^) > 0, 

(3.19) Y{D) < ^.(m)!)"/^^™) _ ^.^^^^^^ 

for large D and using the fact that D — t- Y{D) is continuous, we conclude that 

(3.20) Kn{m):= sup {c„(m)L>"/(2'n) _ y^^^l ^ ^ 

ne[i,oo) 

and 

(3.21) Y{D) < c„(?n)Z)"/(2™) - K,,{m) 

already for all D > 1. Thus, according to ()2.1ip . we have established the following result of [2]. 
Theorem 3.6. The following inequality holds for all u G ^'"(T").- 

(3-22) lkllL(T") < c,,(m)||u||J^'^"^/„"J||(-A)™/2u||2^5„) - Kn{m)\\u\\l2(^jn), 

where Cn{m) and Kn{m) < kn{m) are sharp and are determined by (|3.18p and (|3.20p respectively. 

We see that the constant c„(m) here coincides with the best constant in the analogous in- 
equality on the whole M", see (j3.2p . In contrast to that, the sharp constant Kn{m) is unlikely 
to be expressed analytically for all n and m since, according to the definition ()3.20p . we need 
to know the function V{D) not only for large D, but for all D > 1 in order to compute it. As 
shown in [2] , the asymptotic expansion ()3.17p works for very large D only (if m is large enough) 
and for the intermediate values of D this function is oscillatory (which can be explained by 
studying the limit m — 00, see [2] for the details). 

Moreover, there is a strong difference between the case n = 1 and the multi-dimensional case 
n > 1. In the first shown in [2], we always have Ki{m) > 0, so the lower order term in 
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(j3.22p improves the classical interpolation inequality. However, in the multi-dimensional case, 
this constant become strictly negative for sufficiently large m (for instance, for m > 9 if n = 2 
and for m > 6 for n = 3). Therefore, in that case, the lower order corrector becomes positive 
and necessary for the validity of the interpolation inequality. In other words, if we are interested 
only in the classical interpolation inequality on the torus (without the lower order correctors), 
we have to increase the constant in comparison with the case of M"'. 

Nevertheless, as shown in [2] with the help of numerics there are 3 particular cases where 
Knini) = kn{m) and as a consequence, all constants in (I3.22|) can be found analytically. That 
are 

1) The case n = 1 and m = 1: the inequality 

(3-23) ll^llioo(Ti) ^ ll^llL2(iri)ll^'llL2(Ti) - -|KIIl2(ti) 

holds for all 27r-periodic functions with zero mean; 

2) The case n = 1 and m = 2: the inequality 

^3 241 < ^|U,||3/2 ||,/'||V2 _A||,,||2 

K'^-'^V II^IILoo(TI) ^ 4/2yll^llL2(Ti)ll^ IIl2(T1) 37^ll"lli2(Ti) 

holds for all 27r-periodic functions with zero mean; 

3) The case n = 2 and m = 2: the inequality 

(3-25) II^IIL(T2) < ^ll^llL2(T2)l|Au||i2(T2) - ^lkllL(T2) 

holds for all 27r x 27r-periodic functions with zero mean. As the numerics suggests, that are the 
only cases (at least with integer m) when the constant Kn{m) = kn{m), see [2]. 

Based on the technique developed above, we give below the purely analytic proof of the first 
two inequalities. 



Theorem 3.7. The interpolation inequalities (I3.23P and (I3.24p hold for all 2tt -periodic functions 
with zero mean and all constants in those inequalities are sharp. 

Proof We first consider ^M)- Since $a(x) = (27r)-i/2jr-i((^2 ^ ^^-i^ ^ _i_e-VA|x| aU 
A > 0, according to ()3.14p and summing the geometric progression, we have 

1 7r^/Acoth(7r\/A) - 1 



(3.26) G(A) = - ^ 



27rA 2tt a 

fcez 



Thus, in view of Theorem 12.81 for D > 1 



(3.27) Y{D) = min{(A + D)G{X)} < {(A + D)G{X)}\^^n-i/2 = C^D - 1/2)G{D - 1/2), 

where we have replaced minimum with respect to A by the value at A = Z) — 1/2. Thus, we only 
need to prove that 

(3.28) {2D-l/2)G{D-l/2)-D^/'^ + - = 

= — (iirDV^D - 2 ■ coth a - ■k\/AD - 2 ■ coth a -2- MD^/"^ + A■n^/D] < 

A'k{2D - 1) V / 

for all D > 1, where we set for brevity a := IL^HHEl^ To simplify the expression on the right-hand 
side of ()3.19p . we use that 

y/AD - 2 < 2\fD - -D'^/'^ - — D"^/^ D>1 
2 16 ' - 
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(this inequality is obtained by expanding y^l — {2D)~^ in Taylor series and noting that all 
dropped out terms there are negative). Using also that cotha > 1 in the negative terms, we end 
up with 

(3.29) Y{D) - D^/"^ + -< -r-r-k T (iTTD{8D^^^ - 2D'^^^ - D^^^^/A) ■ cotha- 

vr 4tt[2D — 1) V 

ttVAD - 2 • coth a - 2 - SttD^/^ + AttVd 
(327rL>2(cotha - 1) 



167rL>i/2(2L> - 1) 
8ttD coth a-TT coth a - A^T^/4D - 2\/^ - 8^/D + WirDj < 

< Tii r / 647rZ)_ _ _ ^ _ 4^^4^ _ _ g/p + IGttd] = 

167rL>V2(2L)- 1) Ve2° - 1 7 

~ 167rL>V2(2Z) - 1) l^exp(7rV42? - 2) - 1 ^ 1 + (2D)-i ~ ~ /' 

We note that the function f(D) := — ^ — is strictly decreasing when D > 1 since 

^ ^ exp(7rV4Z)-2)-l o _ 

1287rL> ( - 2 - V4£' - 2 exp(7rV4L> - 2) + vrD exp(7rV4Z) - 2)) 
~ (-1 + exp(7rV4Z) - 2)) V4I) - 2 ^ ° 

for all D > 1. Analogously, the second term in the right-hand side of ()3.29p and the last one 
are also strictly decreasing, so replacing them by the their maximal values at D = 1, we finally 
have 

. s 1/9 1 1 / 647r Stt 

Y(D) - D^l^ + - < — ^ + ^ - vr - 8 

^ ' vr 167rDi/2(2D- 1) Vexp(7r^/2) - 1 2 + ^2 

(-1.3873...) < 0. 



167rL»i/2(2L> - 1) 

Thus, inequality (j3.23p is proved. The fact that all constants there are sharp follows from the 
previously established asymptotics (j3.17p : 

V(D) = dV2 _ i + 0(Z)-i/2) as B ^00. 
vr 

Before we turn to (j3.24p we note that the general formula (j2.44p gives in our case A(-D) = 
D — ^^/D + . . . . Other choices of A in the substitution in (I3.27P will do, for example, A = D — 1, 
X = D — -\^D. However, the choice X = D will not. 

We now consider (13.24p . Using contour integration (see [24l Section 3.3]) we can sum the 
series in G{X): 

1 iry/2 sinh(7r-v/2/x) + sin(7r\/2/i) 1 

^//T^ ~ V" cosh(7r^/2/x) - cos(^^/2^) ~ V' 

and hence 

(3 30) G(X) = A V ^ =]l(^ sinh(7rV2AV^) + sin(^^/2AV^) _ l\ 

^' ' 2TT^^n^ + X vr 1^4A3/4 cosh(7r\/2Ai/4) _cos(^^/2AV4) 2X J ' 



Since 



sinha + sina sinha + 1 -a n AT\i/4 ^ /?r aaa^o 

— < < 1 + 4.1e ° for a = 7rV2A^/* > 7rv2 = 4.4428 . . . 

cosh a — cos a cosh a — 1 
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it follows that 



0(A)<0,A,:.i(^(l.«.-.v...)-±) 

By Theorem 12. 8| for D > 1, substituting A = 3D — 3/2 we obtain 
(3.31) Y{D) = min {(A + D)G{X)} < 

— 1 

< {(A + D)GiX)}\x=3(D-i/2) = i^D - 3/2)G(3D - 3/2) < (iD - 3/2)Gi{3D - 3/2). 
Therefore it suffices to show that 

RiD) := {AD - 3/2)Gi(3D - 3/2) - -^D^/^ + ^<0. 

V 27 37r 

Setting x:={3D- 3/2) ^Z^, x > (3/2) = 1.1066 ... we have 

67rR{D) = {8x' + 3) (l + 4.16"'^^^) - - 2^/27r(x^ + 3/2)1/^ + 4 = 

^ 2.V2 ((. + ^ - ix^ + 3/2)V^) + (. + _L) 4.1e-^^) - ^ ^: ..(.). 

The first term in parenthesis is 0(x~^), hence R{x) < already for, say, x > 2. Instead of 
analyzing R{x) near x = (3/2)^/^ we show in Fig. [T]the graph of R{x) in this region, so that 
R{x) < for all x. 



Figure 1. Graph of the function R{x) on x e [(3/2)^/^^,5]. 

In conclusion we point out that the general formula (j2.44p gives X{D) = 3D-^D^/^ + ..., 

which explains the choice of at least the leading term in the substitution A = 3-D — 3/2 that has 
been used above. The proof of the theorem is complete. □ 

3.2.1.2. The case / > 0. The analysis of the Green's function (j3.12|) for general / > is more 
delicate. Indeed, in this case, the function k — )• ^^2m^x\k\^' ^ singularity at A; = 0. By 
this reason, the associated fundamental solution ^x{x) is not rapidly decaying and the Poisson 
formula ()3.14p becomes not essentially helpful. To overcome this difficulty, following [2], we 
introduce fj, = and rewrite (j3.14p as follows 

(3.32) G(A) = (2')-V E t^,(,|,.p,Lo ^ 1) = "Gl")- 
Then, differentiating the function G{fi) s-times in fx, we get 

(3.33) — ^ = (-l)^r(. + l)(2vr)-" 
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and the function k — )• ^^|^|2(m-i)-|-i)s+i becomes regular at A; = if s > (and the rate of decay 

of this function as A; — >■ oo remains \k\'^"^ for all s). Thus, we may apply the Poisson summation 
formula in order to find the asymptotic behavior of (|3.33p as /i — t- 0+ (analogously to (|3.14p ) 
which after the s-times integration will give us the asymptotics for the initial function G(A) up 
to s integration constants which should be determined using the alternative methods, see below. 

We illustrate this method for the particular case I = 1 and m = 2 only. We also restrict 
ourselves to considering only the two {n = 2) and three (n = 3) dimensional cases although the 
general case can be analyzed in a similar way. 

Let n = 2. Then, we need to differentiate G(/u) only once in order to remove the singularity 
at /e = 0: 

(3 34) dG{^^^^_^^ I 

d/x 47r2/-Jl + /.|fc|2)2- 

Applying the Poisson summation formula to this series, we get 



(3.35) ^ = I ^7 - 1 + 2-V-^/^ E m,{2.^,~'l^\k\) 



where Kn{z) is the modified Bessel function of order n, see [26] for the details. Integrating 
these series in /i and using the standard formulas for the integrals of Bessel functions, namely, 
Kq{x) = —Ki(x), we end up with the desired formula 

(3.36) G{fi) = ^(^logi + ^ + ^-87rE /^o(2^^"^/^|/£ 

where (3 is an integration constant. As shown in [2], based on the Hardy lattice formula 

EnT2i = 4C(i + ^)/3(i + ^), 



where C,{s) and /3(s) are the Riemann zeta and Dirichlet beta functions respectively, 

(3.37) /3 = 7r7 + 4/3'(l) = vr (27 + 21og2 + 31og7r - 41ogr(l/4)) , 

where 7 and T{z) are the Euler-Mascheroni constant and Euler gamma function respectively. 
Thus, we have proved the following result. 

Lemma 3.8. Let n = m = 2 and I = \. Then, at least for A > 0, the Green's function ()3.32p 
can he written as follows: 

(3.38) G(A) = -La- log a + + - ^A- ^ ^o(2vrAV2|^|) < 



k&Zl 



^I^aI^^^s' + ^ + I 

where /3 is defined in ()3.37p . Moreover, the term containing the sum of Bessel functions in the 
middle part of (I3.38P is exponentially small (of order 0{e~'^^^ ^) as A —t- 00. 

Note that the choice n = m = 2 and / = 1 corresponds to the limit case in (j2.9p where the 
extra logarithmic term appears and which is usually called Brezis-Gallouet inequality (see [1]). 
Although the condition (|2.8p is formally violated, most part of the theory developed above works 
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in that limit case as well. In particular, as shown in based on (|3.38p . there exists a constant 
L > such that 

(3.39) Y{D) < ^ {log D + log{l + log D) + L), D > 1, 

47r 

where the constant ^ is sharp and the numerically found value of L = 2.15627. This leads to 
the improved Brezis-Galloet inequality in the form 

(3-40) ||n||^(,.) < ^||V^||i.(,.) (log !!^^!!f^"'^ +log fl + log +l) . 

Let us now consider the 3D case n = 3 with m = 2 and 1 = 1. Then, 

^ ' 87r3 ^ |/fc|2(l + ^|fc|2)' d^ji 87r3 ^ (l + ;u|A:|2)2 

and, applying the Poisson summation formula (|3.13p to the second sum, after the computation 
of the Fourier transform, namely, using (see |22j ) 

(3.42) ^(1/(1 + ^2)2) (^)^_^^H5| 
we end up with 

(3.43) ^ = A - 4^ y e-2.,-/.|.|i/. 
and integrating this series in //, we arrive at 



-l/2g-2^M-l/^|fc|^/^ 



(3.44) G{^i) = ^ j 27rV~'/' + /Sg + ^ - vr 

where (3^ is an integration constant which can be found numerically: 

/33 = -8.91363291758515127, 

see the next section for more details on how to compute it. Thus, we have proved the following 
lemma. 

Lemma 3.9. Let n = 3, m = 2 and 1 = 1. Then, at least for A > 0, the Green's function (I3.32p 
can be written as follows: 

1 / „-27rAl/2|fc|l/2\ 

(3.45) 0(A) = ^ 2,^A'/^ + A + A- - . < 

where f3 is defined in (j3.37p . Moreover, the term containing the sum in the middle part of ()3.45p 
is exponentially small (of order 0{e~'^^^ ^) as A — )• ooj. 

Thus, the Green's function G{X) satisfies the assumptions of Proposition 12.61 with 6 = 1/2 
and gi,g2,g3 = (87r3)~i(27r^, /Sa, 1), respectively. Therefore, 

(3.46) YiD) = ^( in^D^/^ + 2/33 - ~ t""' D^^/A + O(Z)-i), 
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where the third term is negative, which suggests the following inequahty: 
(3-47) ll^xll^ < ^||Vn||i^(Tr3)||A^z||i,(T3) - ^||Vn||i^(Tp3), 

where all constants are sharp. However, up to the moment, we have checked this inequality only 
for large D. The next lemma shows that it holds for all D > 1. 

Lemma 3.10. The inequality ()3.47p holds for every u G -ff^(T'^) with zero mean. 

Proof. We only need to check that the inequality 

(3.48) N(D) < ^ (^47r2Di/2 + 2/3. 

holds for all D > 1. To this end, we will again use the variational representation ()3.27p where 
we put A = D — i and inequality (|3.45|) . Singling out the term 47r^D^/^ + 2/33 this gives 

(3.49) mMD) < (l + jj^) {2.^Vd^2 + /33 + = 

1 / 2 



Att'd'/' + 2/33 + TTTT-^ v^^-f + /33 + 4 + 



2D-iy ^AD{AD - 2) +AD-1 2D - I j ' 

and we only need to check that the last term in the right-hand side is always negative. Indeed, 



2 V4D - 2 2 V4I? - 2 2 1 ^ 

TT — , < TT — , = TT — ; , < 



2 



y/AD{AD - 2) +AD-1 y'AD{AD - 2) +AD-2 VAD + ^/AD - 2 2 + y/2 
if D > 1 and, analogously, 2D-i — ^' Then 

+ ^ + 6 = -0.022892 < 

2 + ^/2 

and the lemma is proved. □ 
3.2.1.3. Computing the integration constants. 

As we have seen above, in the case / > the Poisson summation formula allows to find the 
asymptotic expansions of the Green's function G{X) only up to some integration constants and 
the direct computation of that constants is a non-trivial task since the series p.l2p converge not 
sufficiently fast, especially for big A. Thus, it looks reasonable to find better (e.g., exponentially) 
convergent series for computing them. In the present section, we give an explicit formula for the 
sum (I3.12P in the particular case n = 3, m = 2,l = l considered above in terms of the integrals 
of the so-called Jacobi theta functions and using the known relations for the theta functions, we 
find the formula for the integration constant (3^ through the very fast convergent and convenient 
for computations series. Note also that, although we restrict ourselves to consider only that case, 
the presented method has a general nature and is applicable for computing other integration 
constants including the case of anisotropic tori, etc. 

We consider the Jacobi theta function (see, for instance, [9]) 

(3.50) e3{q)=J2l''- 

Then, the following identity holds which is crucial in what follows: 

(3.51) 03(e-"*) = -^-^3(e-"*"), 

and which follows from the Poisson summation formula ()3.13p with n = 1, f{x) = e~^^/^, 
m = e-«'/2, and = w^. 
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Using also the obvious relation 

1 roo 

(1 - e-*/'^)e-^'* dt, 



1 



k{l + fik) 
we transform (|3.4ip as follows 

f'OO fOO 

STT^Gifi) = (1 - e^*/'^) V [e^'fi[e^'f^[e~'f^ dt = (1 - e"*/^) ([^3(e~*)]^ - l) dt. 
Splitting the interval of integration M+ = [0, 1] U [1, oo) and using (I3.5ip . we arrive at 

/oo 
(l_e-*/M)([03(e-*)]3-l) dt+ 

+ j\l - e-*/^) (7r3/2t-3/203(e-'*-^)3 _ i) dt 

3 /'OO , ^ 

(1 - e-*/^) {Me-')f -l)dt + 7r3/2 y (1 - e-i/^*^))*-!/^ (e^{e-^''f - l) (it- 



Jo Jo 



■^/n dt. 



Now it is not difficult to find the asymptotic expansions for G{n) as — >• 0+. Indeed, as 
elementary calculations show, 

(3.53) / - e-'/^ dt = / i-3/2(i _ g-t) _ / ^-3/2(;^ _ g-t) ^ 



-1/2/ 



/"OO 

(2V^ - / r3/2 + 0(e-i/^)) = 2V^/x-i/2 _ 2 + 0^(1). 



Using now that O^i^x) — 1 = 0{x) as x — )■ we can pass to the limit /x — )• in the integrals 
in ()3.52p containing 63 and comparing with p.44p . we see that 

/oo pCO . , 

(^3(e-*)3 _ 1) dt + ^3/2 i-1/2 (^^3(e-^*)3 _ l) dt. 

We do not know whether or not the integral in (j3.54p can be computed in closed from, however, 
it is convenient for high precision numerical computation of the constant /33. Indeed, expanding 
back the Jacobi function in Taylor series and using that 

Jx Vk 
where erfc(x) := e~*^ dt, we end up with 

(3.55) ft = -l-2,3'^+Ef^? + -*"""^ 



|A;|2 |A:| 



We see that the rate of convergence of the series is super-exponential and using Maple to compute 
it, we get the desired value 

(33 = -8.91363291758515127 
which has been used in the previous section. 
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3.2.2. Inequalities on spheres. We recall the basic facts concerning the spectrum of the Laplace- 
Beltrami operator on the {d — l)-dimensional sphere S*^"^: 

-AY^ = AnY^, k = l,...,kd{n), n = l,2,.... 

Here the are the orthonormal spherical harmonics. Each eigenvalue 

An = n{n + d-2) 

has multiplicity 

, , , 2n + d-2fn + d-3 
kd[n) = ■ 



n \ n — 1 

In particular, for d = 3, 4 we have 

: An = n(n + 1), /csfn) = 2n + 1, 

(3.56) 

§3 : A„ = n(n + 2), k^n) = {n + 1)^. 
The following identity is essential [22|: for any ^ G S"^"^ 

(3.57) Y^iif = 

where a{d) = 2tt'^/'^ /V{d/2) is the surface area of S'^~^ . 

Finally, since the kernel of A is the one-dimensional subspace of constants, throughout below 
we assume orthogonality to constants: 

H'{S^-^) = {^€ H'{S'^-^), [if, 1) = 0}. 

3.2.2.1. Inequalities on spheres: S^. We consider below applications of the general theory to 
inequalities on the 2D sphere with A = (—A)*", m > 1/2, and B = I, where A is the Laplace- 
Beltrami operator. The first positive eigenvalue of —A is 2 and in accordance with ()2.13p we 
consider the operator 

A(A) = (-A)'" + A/, A>-2'". 

Its Green's function Gx{x,^) is 

n=l k=i ^ ^ " 

and thanks to (j3.57p the function Ga(C>C) is independent of ^ and is given by 

00 2ri+l Vklt\l 1 °° 9r( -U 1 

(3.58) G(A) = GA(e, i) = yy , , r = , , r . 

^ ' ^' ^ ^ ^1 + 1))™ + ^ 47r (n(n + 1))- + A 

Setting n := A~^/™ we have 

-. 00 

G(A) = y(.2n + l)<^(/.n(n + 1)), 

where 



4tt 

n=l 



V9(rc) 



x"^ + 1 

Since ^{0) = 1 and f'{0) = 0, Lemma 13.141 below gives 



(3.59) G(A) = i- 



A-^i^ - ^A-i + • A-(2-9) 
3 



+ 0(A 



-(3-26»)^ 
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where 6 = (m — l)/m and 

K = / ip{x)dx = — — . 

In other words, we have shown that for /g(A) = G(A) the asymptotic expansion (|2.4ip holds 
with the fohowing g2, gz independent of ^ € S^: 

K 1-e 1 

(3-60 g^ = — = , 52 = -7-1 53 = 0. 

Att 4 sm an dtt 

Next, differentiating (j3.58p and again using Lemma l3.14l we obtain the asymptotic expansion 
ior g^{X) = g{X): 

g{x) = -/'(A) = -G'(A) = — y ,^ , ^''^^ = —ij?^ y , = 

' ^' 47r ^ ((n(n+ l))'" + A)2 47r^ ((;un(n + 1))™ + 1)2 



1 

47r 

where we used 



A-(i-^)ejf-^A-2 + 0-A-(3- 



+ 0(A 



-{4-e)^ 



^ ' Jo + 

This justifies differentiation of the asymptotic formula for /(A), as required in Proposition 12.61 
Applying it we obtain as a result the asymptotics of the function Y{D) on S^, 

Theorem 3.11. The function Y{D) solving on tfi^ extremal problem 



(3.61) 



V(^,D) :=sup||u(OP : ueH"\S^), \\uf = 1, ||(-A)'"/\f = Z)| 



is independent of § and has the following asymptotic behavior as D ^ oo: 
where 9 = (m — l)/m. 

We now apply the asymptotics ofY{D) obtained above to multiplicative inequalities with re- 
mainder terms on S^, gy Theorem 1 2 . 5 1 the sharp constant K = Kg in the classical multiplicative 
inequality 

(3.63) u{^f < \\u\\l^ < A'0||uf^||(-A)'"'uf (1-^) 
is given by 

(3.64) ^« = «r^S^"°<^' = s^T^Si;'"£ 



(n(n + l))'" + A' 



We note that (I3.64p was obtained in [12] by a somewhat similar but less general argument than 
the one used in Theorem 12.51 It was also shown there that for (integer) m, 2 < m < 7, we in 
fact have that the supremum is attained at infinity 

dx 



poo 

(3.65) supA^G(A) = lim A^G(A) = / - 

A>0 Jo 1 

which gives 



+ X'' 



K. ' 



4 sin Ott 
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and, equivalently, 

(3.66) V(z,).v„.(0)<j-i^(i^)^'- 

for 2 < m < 7. However, for larger m's the supremum in (|3.65p is attained at a finite point 
A* < oo. An explanation of this phenomenon for the torus has been given in [5]; below we 
consider the case of the sphere S^. 

Lemma 3.12. For all sufficiently large m the function h{X) := A^G(A) attains a global maximum 
at a finite point A* and /i(A*) > h{oo). 

Proof. Setting in ()3.64p A = z^^™ we see that up to a constant factor, h{X) is equal to 



2n+l 



n=l 

We consider the following partitioning of the half-line x > by the points 

(n — l)n 

On = an[i^) = 5 , n = 1, . . . . 

Then a direct inspection shows that 

ui \ ^ ( \( ^ , ^ (/?(an) + v?(an+i) , ^ 
(3.67) li[y) = -99(a2)(a2 - oi) + 2^ ^ (a„+i - a„), 

n=2 

where (^{x) = 1/(1 + x'"), which looks like a step function for large m: (^(2^) ^ 1 for x < 1 and 
(/^(x) ~ for X > 1. In view of Lemma |3. 141 below we have 

i7(oo) := lim }i{v) = / ip{x) dx = — = 1 + Om^oo(l)- 

u-^oo Jq sinvr/m 

We fix a large m, and set, say, v = vq = \f2 + 1/100. Then 

ai(i^o) = 0, 02(^^0) < 1(= 0.986), 03(1^0) ~ 3(= 2.958); 

02(1^0) - ai(i^o) = 0.986; 03(2^0) - 02(1^0) = 1-972. 

Since •^(0.2) = 1 + Om-5-oo(l)i ^ifl.'i) = Om^oo{^), and the sum from ?i = 3 to 00 in (|3.67|) gives a 
contribution of the order Om,--s.oo(l)) it follows that 

986 1 972 

H{i^o) = + -Y~ + (''n^ooil) = 1-479 + O„^oo(l)- 

Hence H[uq) > H{oo), and the proof is complete. □ 

This argument also explains the initial oscillatory behavior of H^v) when the corresponding 
node anil') hits 1 with the increase of ly. Figure [2] shows the monotone behavior of i^(i^) for 
m = 2 and the initial oscillations of Hii/) for m = 10. Observe that the first (and global) 
maximum in the second case is located near z/ = ^/2. 

Turning to the inequalities with remainder terms we note that the third term in (|3.62p is 
negative, therefore the improved inequality 

(3-68) Y{D) < —1— (^—^) ' D'-' - 77^ V(D) 

holds for all D > Dq, where Dq is sufficiently large. On the finite interval p^jDo] computer 
calculations are reliable; their results are shown in Figl3l 
In particular, for m = 2 

Y{D) < Y{D) for all D, 

which proves the following result. 
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Figure 2. Graphs of the normahzed function H{i') for m = 2 and m = 10 





Figure 3. Graph of the function Y{D) - Y{D) for m = 2, 3, 4. 



Theorem 3.13. Ifu(^H\E>^), then 



(3.69) 



1 



1 



< —\\u\\ IIAull \\u 



loo - 



where both constants are sharp and no extremal functions exist. 

For m = 3 and m = 4 the function Ym{£,, D) —Ymi^,, D) attains a global maximum 0.00486 . . . 
at -D = 15.8 . . . for m = 3 and a global maximum 0.0189 ... at D = 22.4 . . . for m = 4, 
respectively. Accordingly, we have 



(3.70) 



1 



'""""^ - 4sin^^ 



u 



1 26 I 



0™/2n||2(i-e) 



where m = 3, 4, ^ = (m — l)/m and £3 = 0.938 . . . , £4 = 0.821 .... 



We now prove the asymptotic formula for functions defined by the following series (|3.7ip that 
we have been systematically using above. Let F{n) be defined as follows 



(3.71) 



F(^) :=^(2n+l)/(^n(n+l)), 



n=l 
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where / is sufficiently smootli and sufficiently fast decays at infinity. We need to find the 
asymptotics of F{fi) as /i — )• 0. 

Lemma 3.14. The following asymptotic expansion holds as ^ ^ 0; 

1 r°° 2 1 

(3.72) = -j^ f{x)dx - -/(O) - -^/'(O) + 0{^?). 

Proof. We set 



R{x) = (2x + !)/(// + 



and observe that 



(3.73) / R{x)dx = - f{fix{x + l))d{nx{x + l)) = - f{x)dx. 

Jo Jo Jo 

We calculate the derivatives of R up the order 5 at x = 0: 

^(0) = /(O); 

i?'(0) = 2/(0) + ^*/'(0); 

/?"(0) = 6/i/'(0)+^V"(0); 
^^■^^^ ^'"(0) = 12^/(0) + 12^.V"(0) + ^V'"(0); 

i?{4)(o) = 60^V"(0) + 20^V"(0) + //^^^(O); 

i?(5)(0) = 120^V"(0) + 180;uV"(0) + 30//(^)(0) + mV^^^(O), 

and, in addition, 

(3.75) R^^Hx) = 840/i^(2x + l)f^^\fix{x + 1)) + 420/(2x + iff^'^^fixix + 1))+ 

+ 42^5(2x + iff^'^^fxxix + 1)) + f,\2x + lYf^'^Hpxix + 1)). 
Next we use the Euler-Maclaurin formula (see, for instance, |16| ) 

(3.76) Vi?(n)= r R{x)dx + -R{0)-Y^R'''~^\^)- H ^4r^R^''\x)dx, 
fo Jo 2 ^ i\ Jo k\ 

where the B^s are the Bernoulli numbers: B2 = ^,,63 = 0, = — -B5 = 0, i?6 = ^) • • • i 
and the i?fc(x)'s are the periodic Bernoulli polynomials. Using (j3.76p with k = Q and taking into 
account (I3.73P and (I3.74p we obtain 

F{^i) = -R{0) + Rin) = - / f{^)dx - ^/(O) - -(2/(0) + m/(0))+ 

n=0 ^-^0 2 ^2 



+^(12/z/'(0) + 12/.V"(0) + ^^V'"(0))- 
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-L_(120^V"(0) + 180^.3^'"(0) + 30//(')(0) + ^.^('^(0)) - ^R(^){x)dx. 



This gives (j3.72p provided that the remainder integral term is of the order 0{p?). The periodic 
Bernoulli polynomials are clearly bounded on (0, 00). Therefore the contribution of each term 
in (j3.75p is of the order /i^. For example, the last term is of the order 



00 1 roo 



1^1^ I x'^g{fix^)dx = ^fi^ I y^g{y)dy. 
Jo ^ Jo 

The three remaining terms in (|3.75p are treated similarly. □ 
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3.2.2.2. Inequalities on spheres: S'^. We consider on the 3D sphere S'^ only one example with 
/ = 1 and 7n = 2, so that 9 = 1/2. We set A = (-A)^, B = -A, and let 



The Green's function of A(A) is 



A(A) = (-A)2 - AA. 



(n+l)2 



:[ fit n(n + 2)(n(n + 2) + A) 
and again using (I3.57P we see that Ga(C,0 independent of ^: 



(3.77) 



G(A) = GA(e,e) = ^E 



(n+l)2 



for any ^ G S'^. 



27r2 ^ 2)(n(n + 2) + A) 

We do not need an analogue of Lemma 13.141 here since the series in (|3.77|) can be summed 
exphcitly [12] : 



(3.78) 
and, in addition, 
(3.79) 



G(A) 



1 VA-lvr ^ r-^ T 1 

^^2^°'^^^^-^+4A 



supA^G(A) = lim A^G(A) = -3- 
A>o 47r 



which gives in view of Theorem 12.51 the following sharp inequality on jl2j 



(3.80) 



\u\\i, < — ||Vn||||An| 



To obtain power expansion of G{X) we can replace coth vrV A — 1 by 1 which gives : 



= 2^ 



so that (f2:in) holds with: 

(3.81) gi = — , 92 

Using Proposition 12.61 we obtain the following result. 
Theorem 3.15. The solution Y{D) of the extremal problem 



3 

'8^' 



(3.82) 



V(C,Z)) :=sup<' |u(OI^ : u G H'^{S^), ||Vuf = 1, \\Au 



D 



is independent o/ ^ G and has the following asymptotic behavior as D ^ oo: 



2tt 47r2 



(3.83) V{D) 

The third term in (|3.83|) is negative, therefore 
(3.84) 



l + 47r' 
167r3 



Z)-i/2 + 0(D-i). 



for all D > Dq, where Dq is sufficiently large. However, similarly to Theorem l3 . 71 and Lemma [3.101 
taking the advantage of the explicit formula (|3.78p . we have the following result. 

Lemma 3.16. Inequality p.84p holds holds for all D > -y/S. 
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Proof. In view of Theorem ESI for D > y/3 (the first eigenvalue of (-A)^/^ on ig ^) 

(3.85) Y{D)= mm {{X + D)G{X)} < {{X + D)G{X)}\x=D = 'iDG{D). 

Hence, using ([37781) . the inequahty VD-l < /D- 1/{2VD) - 1/{8D'-^/^), and replacing cotha 
in the negative terms below by 1, we obtain 

Y{D) - + < 2DGiD) - ^Vd + = 

2D /y/D-lTT , D-1 1\ \ I— 3 

TTT ^-F. — TTCothWl)- 1 ^ + tt: - T^vl) + = 

i {sfD (co.h ,yirrT-i)--i=-_L_ + -i^) = 

since Fi(L») = (cothvrVi^ - 1 - l) is decreasing for D > ^3 and Fi(L>) < Fi(^/3) = 
3coth(7r\/V3 - 1 - 1) = 0.0278. . . , while F2(L>) = Z)/2 + 1/8 - 2^JD|^: is increasing and 
FiiP) > F2{y/3) = V^/2 + 1/8 - 2^^/tt = 0.153 .... The proof is complete. □ 

As an immediate corollary we obtain inequality ()3.80p in the following refined form. 

Theorem 3.17. If u € H'^{S^), then 

(3.86) \\u\\l < i-||Vn||||An|| - ^llVnf, 
where both constants are sharp and no extremal functions exist. 

Remark 3.18. The coefficient of the leading term in the asymptotic expansions ()3.62p . (I3.83P 
coincides with the corresponding constant in M". 



3.3. Manifolds with boundary. 

3.3.1. First-order inequality on the half-line. We consider a couple of inequalities on the half-line 
IR+ = (0, oo). On the whole M we have 

(3.87) hllL < Ikllh'll, 

where u G ff^(IR) and where the constant 1 is sharp and there exists a unique extremal function 
ti*(x) = e~I^L Using the extension by zero we see that inequality (|3.87p still holds on M+, and 
there are no extremal functions since ti*(x) > 0. We can obtain the following refined form of 
this inequality 

(3.88) u{Cf < ||u||||n'||(l - e-2«), u € H^{R+) 

from the general theory developed above. In fact, the Green's function of the operator 

Au = —u" + u, u(0) = n(oo) = 



IS 



G(x,e) 



e ^ sinhx, x < ^, 
e~^ sinh^, x > ^. 



As in Lemma |2. II we have 
and (I3.88P follows by the standard scaling argument 



< [WuT + \\uf)G{i,i) = [Wut + \\uf)^-^ 
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3.3.2. Inequality for Bessel operator on the half-line. Here we consider a one-dimensional in- 
equality on the half-line with Au = —u" — n/4x^ and Bu = u. Although the corresponding 
sharp constant was found in [T7] in connection with Lieb-Thirring inequalities for radial poten- 
tials (see also [8]) we include this example to illustrate our general approach. We first observe 
that in view of the one-dimensional Hardy inequality 

^^dx < 4 r{u'{x)fdx, 
X Jo 

the operator A is non-negative. In accordance with the general theory developed in section [2] in 
order to find the sharp constant C = C(^) in the inequality 

(3.89) \u{0\' < C{0 (^f [{u'{x)f - ^) dx^ 

we need to write down the Green's function of the following operator 

A{X)u = -u" - u/Ax'^ + Alt 
with Dirichlet boundary conditions n(0) = u{oo) = 0. We have 

fo Qr,^ r VxlKo{V\i)lQ{y/Xx), X < 

""^^"'^^ - i v^cKo{Vxx)io{Vxo, X > e, 

where Kq and Iq are the modified Bessel functions of zeroth order. In fact, the functions 
^/xKo{y/Xx) and ^Io{VXx) both satisfy the homogeneous equation, /o(0) = 0, Kq{oo) = 0, 
and the jump condition G(^ + 0, ^) — G(^ — 0, ^) = 1 is satisfied in view of the Wronski identity 
Iq{x)Kq{x) — Iq{x)Kq{x) = —1/x. Therefore by Theorem 12.51 the constant C(^) is, in fact, 
independent of ^ and is given by 

C{0 = 2sup^/AGA(e,0 = 2snp^^Ko{VxOlo{VXO = 

A>0 A>0 

= 2supri^o(r)/o(r-) = 2r/^o(r)/o(r)|r-=r.=i.075... = 2-0.533-- - = 1.06--- =:a. 

r>0 



Furthermore, for every fixed ^ > inequality ()3.89p turns into equality for C{^) = C* and 
u{x) = Gx^{x,^), where A^, = A*(^) = r^/^^. We also observe that Io{x) ~ 
Ko{x) ~ \/^ as X — )• oo, and hence 



lim rKo{r)Io{r) = ^. 

3.3.3. First-order inequality on interval. In this section we consider the correction term for the 
interpolation inequality 

(3.91) ||it||^ < 1 • u£H^{0,L) 

on an interval (0, L) with zero boundary conditions. The constant 1 is sharp, since it is sharp 
for the inequality (|3.87p on M and we can use extension by zero. The unique extremal function 
in (|3.87p does not vanish, therefore there are no extremal functions in (j3.9ip . and one might 
expect that an inequality similar to ()3.22p holds in the case of a finite interval: 

(3.92) ll-ull^ < ||n||||n'|| - CL||nf, ueH^{0,L). 

In fact, cl = as we now show. First, by scaling, if (j3.92p holds, then cl = c/L where c is 
an absolute constant. Next, we consider the truncated extremal function 

< X < a; 

<Pa{x) = { e~'^(a + 1 — x), a<x<a + l; 
0, X > 0, 
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and set ipa{—x) = ifaix). Then 

II, o II 1 II, o l|2 1 -'■0-21 11,0' l|2 1 I ^-2a 

Substituting this into (|3.92p and letting a — )■ oo we obtain that (j3.92p can only hold with c = 
as claimed. 

However, the correction term exists but is exponentially small. More precisely, we have the 
following result. 

Theorem 3.19. Let u G H^{0,L). Then 

L\\u'\\ 

(3.93) lltill^ < ||n||||n 11(1 - 2e 

The coefficients of the two terms on the right-hand side are sharp and no extremal functions 
exist. 

Proof. Without loss of generality we set L = 1. The Green's function of the boundary value 
problem 

-y" + Xy = 6ix,0, 2/(0) = y(l) = 

is 

(o OA) r F\- 1 / sinhAV2:^sinhAV2(i _ ^), < rr < 

^d.y4j ^^^'''^^ - Ai/2sinhAi/2 \ sinhAV2^sinhAi/2(l-x), ^ < x < 1. 

In addition, G\{^.,^) attains its maximum value with respect to ^ at ^ = 1/2: 

As in Lemma l2. II 

u{if < GxiC,mWf + X\\u'f) < Ga(1/2, l/2)i\\uf + X\\u'f), 
where we have the equality at the point 1/2 for u{x) = const Ga(x, 1/2). Hence, 

Y{^,D)<Y{1/2,D), 

and as in (j2.30p . ()2.3ip we obtain the parametric representation of Y{D) := V(l/2, D): 

(3.95) ^(^) = ^' V{D{X)) = l^, AG[-vr2,oo), 

5(A) g{X) 

where 

(3.96) /(A) = ^ tanh ^(A) = -/'(A), h{X) = /(A) - A<7(A), 
and (|3.93p is equivalent to the inequality 

(3.97) V(D) < Vd{1 - 2e-^) for D > vr^. 
Using the variational representation from Theorem 12.81 we have 

(3.98) V(D)= min J_ tanh ^ • (A + D). 

Ae(-7r^oo) 2VA 2 



To simplify notation we denote d:= yD and further set 

(3.99) A* = X^{D) = D{1 - 2VDe-^f = d\l - 2de-'^f. 

As we shall see below A,,, {D) contains the first terms of the asymptotic expansion as Z? — t- oo of 
the unique solution of the first equation in (j3.95p . or, equivalently, of the unique point A where 
the global minimum in (j3.98p is attained. For the moment, however, we just substitute X^{D) 
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into the right-hand side of (j3.98p and see that we prove (j3.97p if we can show that the following 
inequality holds for all d > vr: 

(3.100) L(l-2&-^) tanh ^^'"f^"'^ {d\l-2de-r + d')<d{l-2e-% 
or 

(3.101) tanh ^(^"^^^ [(1 _ 2de-'^f + A < 2(1 - 2e-'^)(l - 2de-'^). 



2 

Next, we use that tanh(ti/2) < 1 — 2e~" + 2e~^", and we also observe that the quadratic 
polynomial 1 — 2x + 2x^ = 2{x — 1/2)^ + 1/2 is monotone decreasing for x G [0, 1/2] so that 

tanh ^(1-^^^""') < ^_2e~d(i-2de-^)^^^-2dii-2de-^) < i_2e-'' {l+2d^ e-'')+2e-^'' {l+2d\-'')' , 

where we used that e~'^(^~^'^'^ ''^ > e~'^{l + 2d'^e~'^) which essentially is the inequality > 1 + x. 
Combining the above we see that it suffices to establish for d > ir the inequality 

F{d) := (l - 2e-'^{l + 2d^e-'^) + 2e-'^\l + 2d^e-'^f^ [il - 2de-^f + l) - 

-2{l - 2e-'^){l - 2de~'^) < 0. 

Simplifying F{d) we obtain 

F{d) = 4e-2'^Fi((i), 

where 

Fi{d) = -d^ + 2d^e-'^ + U^e-'^ + 1 - 2de-'^ + 2d^e-'^'^ - 8d^e-^'^+ 

+8d^e-^'^ - M^e-^'^ + Sd^e"^"^. 
Next, dropping all negative terms except for — we have 

Fi{d) <-d^ + 1 + F2{d), F2{d) := 2d^e-'^ + Ad^e''^ + 2^^6-2^ ^ §^4^-3^ ^ §^6^-4^^ 
Each term in F2{d) is monotonely decreasing for d > ir, and hence 

Fi{d) < -7r2 + l + F2(7r) = -2.530... < 0, 

which proves ()3.97p . 

To explain our choice of A=i,(D) in ()3.99p we find the asymptotics as — >• oo of the inverse 
function A = X{D) since our A*(L') contains the first two terms of this asymptotic expansion. 
We have 

/(A) = ' tanh ^ = J^',^^ = W^l^ (l - 2e-- + 2e--^^^ + 



2A1/2 2 2A1/2 1 + e-^'/' 2 

g(A) = iA-3/2 (l + (-2 - 2Ai/2)e-^''' + (2 + A\^l^)e~^'^''' + 
h{\) = 1a"1/2 (l + (-2 + 2A^/2)g-AV2 ^ ^2 _ 4Ai/2)e-2^^/' + 
and 

^ MA) ^ A(l + (-2 + 2AV^)e-^^^' + (2 - 4AV^)e-2^^^' + . . . ) ^ 
^ ^ 5(A) i + (_2-2AV2)e-Ai/2 + (2 + 4AV2)e-2AV2 + , ^ 

= A + 4A3/V^^/'+8A2e-2A^'' + .... 
Hence the inverse function \{D) has the asymptotics 

\{D) = D{1 - 4/De-^ + ...) = d^{l - 2de-'^ + ■■■?, 
whose first two terms give (j3.99p . 
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Finally, 

YiD) = — tanh ^ ' ■ {X{D) + D) 

1 1 , d(l-2de~'^ + ...) / ,2, , ,2 ,2 

2 d(l - 2fie-'^ + . . . ) 2 V ^ ^ 

d.tanh^fc^^^|I^-±^(l-2de-'^ + ...) (l + 2cie-^ + . . . 

d. tanh fc^^|I^^±^(l + ...) = d(l-2e-'^ + ...), 

which proves sharpness and completes the proof. □ 

Remark 3.20. Arguing similarly to Proposition 12.61 one can write down the three-term expan- 
sion of Y{D) as D ^ oo: 

(3.103) YiD) = D^l^ - 2D^l^e~^'"" - 2D^I\D - l)e"2^"''' + o{e~^''-"'+'^). 

The third term is negative and hence inequality (I3.97P holds for all sufficiently large D. However, 
as we have shown, (|3.97p holds for all D > ir^. 

3.3.4. On a second order inequality on the interval. We want to apply the above developed 
theory to the following one-dimensional interpolation inequality 

(3.104) ||n||^ < u G ^^(0, 1) n ^^(0, 1) 

Here, Au{x) := u^'^^x) with boundary conditions u{0) = u{l) = u"{0) = u"{l) = 0, B = Id, 
9 = 3/4 and A(A) = u^'^^ + \u and, in order to find the key function G\{^,^), we need to solve 
the equation 

u""{x) + \u{x) = 5{x-i). 
Using the orthonormal system of eigenfunctions {-v/2 sinvrnx}^]^ we obtain 

oo . . ^ 

^ , , v-^ sm vrna; sm vrnt 

(3.105) GA(x,g) = 2^ ^^^^^^ S 

71=1 

and setting A = ^o^tt^ to simplify the formulas below we have 

^ oo . 2 /■ 

n (c c\ ^ sm^7r< 
^A(C,e) = ^2^;^T^. 

n=l 

Next, the identity sin^vrn.^ = 1/2 — (e^'^*"? + e~^'^*"^)/4 and the Poisson summation formula 

oo oo 



give 



sin^vme _ 1 / ^ 1 ^ e^-'< 



n=l 
1 1 



47J^[ E /(^^«27rn)- ^ f{V2a27r{n + 0) 

' \n=— oo n=— oo i 

1 1 7r^/2 



^^n=— oo n=— oo 

oo oo 



4(V2a)3 2 



Y /o(a27rn)- ^ /o(«27r(n + 0) 
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where 



fiy) 



fo{^), fo{x) = e '""'(cos |x| + sin|x 



Combining the above, we obtain for the key function in Theorem 



(3.106) X'/^Gx{^,0 



V2 



J2 Ua27Tn)- J2 /o(a2^(n + 0) 



V2 



We now observe that the function fo{x) (which up to a constant factor is the fundamental 
solution of the operator (I3.3.4p on the whole line) has positive local maximums at x = 27rn (the 
one at being the global) and negative local minimums at x = vr + 27rn (the ones at zLtt being 
the global). 

For a large enough and ^ < 1/2 the leading terms are 

Sia,0 = MO) - /o(2vr<) + ©(e'^-'^), 

and since fo{x) has a negative global minimum at vr, we see a "horn" of height at ^ = a~^/2, 
which gives the maximum value of S{a, ^) for large a; see Fig. HI 

As for the global maximum of S{a,^), we see that if a = 1 and ^ = 1/2, then the first and 
the second sums in (I3.106P count one by one all the maximums and all the negative minimums 
of /o, respectively. Therefore 

S{a,0 < 5(1,1/2) = coth|. 
Therefore we have proved the following result. 
Theorem 3.21. The sharp constant in inequality (j3.104p is 

V2 



K 



)^ • coth — 
^ 2 



1.09033... . 



The unique extremal function is given by (|3.105p with X = 4tt'^ andi = \. 

Remark 3.22. We point out that the sharp constant in (j3.104p for u € Hq{0, 1) is the same as 
that on the whole line, namely, 

V 27 




Figure 4. Plots of ^ S{a,^) for a = 2.5/7r (red), a = I (green), a = A/tt 
(brown), a = S/vr (blue), a = (7, 15,20)/7r 
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3.3.5. A second-order inequality in 3D. In conclusion we consider a second-order inequality 
in a three dimensional domain $7 C £qj. -v^hich the passage from Hq{Q) to a wider space 
H'^{Q) n Hq{Q) does not increase the constant in the corresponding interpolation inequality. 
This inequality was obtained in [2^ by a somewhat problem specific method, nevertheless we 
present its proof in the framework of our general Theorem 12. 5[ 



Theorem 3.23. Let C arbitrary domain. Let Hq{^) be the completion of C^{^) 

in the norm ||V • HljCH)- Then for u € Hq{^) H {u : Au £ L2{^)} the following inequality holds 

(3.107) II^IlL < ^||Vn||||An||, 

where the constant is sharp and no extremal functions exist, unless Q = W^ and 

uix) = T-\i/m' + \^\')). 

Proof. We first assume that is a bounded domain with smooth boundary. Then by the elliptic 
regularity we have 

u G H^{n) n{u: Aue L2{n)} = H^{n) n H^{n) 

In accordance with Theorem 12.51 we have to consider the Green's function Gx{x,(^) of the 
following 4th order elliptic equation 

.3 ^Qg. (A' - AA)GA(x,e) = -A(-A + X)Gx{x,0 = S{x - 0, 

Gx{x,0 = AGx{x,0 = 0, C (^fl, X £ dn. 
We denote by gD{x,£,) the Green's function of the Dirichlet Laplacian: 

-AgD{x,0 = S{x-0, gD{x,O = 0, C^^, xedn, 
and by gH{x,Cj the Green's function of the Helmholtz equation: 

(-A + \)gH{x, i) = 5{x - 0, gnix, = 0, C £ x € dn. 
By the maximum principle we have 

(3.109) o<_g^(^x,0 <—. TT, 0<gH{x,O<-ri ^' 

where the functions on the right-hand sides are the corresponding fundamental solutions in 
Therefore 

f f e~^l^'~^l 1 

0<Ga(x,O= / gH{x,y)gD{y,£,)dy < / r—. Tjdy < 

Jn Jq 47r|x - y\ 47r|y - ^| 

f e-^l^-S'l 1 1 f e*(«-^>^(iz 

< / ^ — 7^(^y 



47r|x - y| 47r|y - ^1 (27r)3 Jiga |z|2(|z|2 + A) ' 

the function on the right-hand side being the fundamental solution of ()3.108p in W^, and, con- 
sequently, 

. . _J_ f dz ^ 1 

^Al^,U< (2vr)3yj,3 |z|2(|z|2 + A) 47rVA' 
Now (j3.107p follows from Theorem 12. 5i The constant l/(27r) is clearly sharp, since we have the 
inequality with the same constant for u G Hq{^)^ see ()3.2p . and (using extension by zero) for 
u G Hl{Q) C Hl{Q) r\{u: Au £ L2{n)} . 

Exhausting an arbitrary Q C with bounded domains J7„ one can show that (j3.107p holds 
in the general case (see [27] for the details). □ 

Remark 3.24. Unlike the previous example, the fundamental solution of (j3.108p in 
positive. 
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